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ABSTRACT 

The  approximation  of  analytic  or  harmonic  functions  on  plane  finitely 
connected  domains  with  smooth  boundary  curves  is  considered,  and  several 
easily  computable  families  of  near-orthonormal  functions  are  described.  The 
functions  can  be  used  in  fast  Galerkin  schemes  for  numerical  approximation. 
Several  numerical  examples  are  presented. 
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functions,  approximation. 
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SIGNIFICANCE  AND  EXPLANATION 


In  many  applications  one  Is  Interested  in  approximating  analytic  or 
harmonic  functions  on  finitely  connected  regions.  The  families  of  analytic 
and  of  harmonic  functions  described  in  this  paper  have  the  following  features 
in  common  which  make  them  attractive  for  efficient  numerical  computation: 

1 )  They  are  easy  to  compute . 

2)  They  are  near-orthonormal  with  respect  to  an  inner-product  on  the  boundary 
of  the  region. 

3}  They  permit  rapid  approximation  schemes.  Only  0(n  log(n) )  operations 
are  required  as  n  where  n  is  the  number  of  independent  functions 

in  the  approximant. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  HEC,  and  not  with  the  author  of  this  report. 


2.  Near*orthonormal  functions  on  regions  with  bounded  simply  connected  complement » 


We  present  three  families  of  near-orthonormal  analytic  functions  for  the  case  when 

n  is  the  complement  of  a  simply  connected  bounded  closed  region  P  ;=*  t\5l  with  smooth 

c  ^ 

boundary  3Q.  (t  denotes  the  extended  complex  plane.  The  functions  defining  these 
families  are  of  the  form 


(2.1) 


q^(z) 


Pk<^’ 


I 

r 

1=1 


k/t 

( z-w , ) 

1 


k  =  0,1,2,... 


'  Pk'*’ 

(2.2)  q^(z)  :=  '-J -  ,  k  =  0,1,2,... 

n  ( z-w . ) 

.=1  3 


(2.3) 


qj^(z)  :=  rational  function  wit)i  all  poles  In  (1^,  k  =  0,1,2,...  , 


w)iere  P|,(z)  denotes  a  polynomial  of  degree  <  k,  and  t)ie  w^ts,  j  =  1(1)t,  are  points 

1/1 

in  Q^\30,  not  necessarily  distinct.  Tlie  branch  of  the  root  (z-w^)  is  selected  so 
that  its  argument  varies  continuously  as  z  traverses  90.  The  w^is  are  defined  by 


*'y  1  *  1(1)*,  k  =  0,1,2,... 


To  describe  the  function  feunilies  more 


the 


precisely,  first  select  points  w.,  j  =  1(1)t,  in  $J  \9f*.  Let  G(z,w.)  denote 

J  c  J 

Green's  function  for  the  Laplace  operator  for  with  a  logarithmic  singularity  at  w^, 

i.e.,  &G(z,w,)  »  0  for  z  €  0  \{w,),  G(z,w^)  =  0  for  z  e  3S1,  and 

3  c  1  3 

I  •:r~  (z,w.)|dzl  =  1,  where  -s—  denotes  the  normal  derivative  into  U  .  Then 

3(3  3"  3  3n 


(2.4) 


1  *  3r 

J(z)  !=  7  ).  ■JT  (Z.W.).  z  e  30  . 
*  3n  J 


Let  C  be  an  arbitrary  but  fixed  point  on  30,  and  define  for  z  B  30, 
o 


SOME  FAMILIES  OF  NEAR-ORTHONOBMAL  ANALVTIC  AND  HARMONIC  FUNCTIONS 


Lothar  Reichel 


1 .  Introduction 

Let  n  be  a  finitely  connected  open  region  in  the  extended  complex  plane  with 
boundary  which  consists  of  smooth  Jordan  curves*  We  consider  the  approximation  of 

functions  f(z),  analytic  in  !1,  smooth  on  0;  =  f!  U  3!J,  and  known  on  3D.  By 
potential  theoretic  methods,  we  derive  several  easily  computable  families  of  functions 
which  are  analytic  on  D,  and  are  near-orthonormal  with  respect  to  an 
inner-product  defined  on  SD.  Also  the  functions  {Re  Re  q^(z),  Im  q^(z). 

Re  q2(z),  im  q^fz),...}  are  near-orthonormal  on  3D,  and  therefore  suitable  for 
approximation  of  harmonic  functions.  The  approximants  are  computed  by  a  Galerkin 
scheme.  After  some  Initial  computations,  which  do  not  have  to  be  repeated  for  each 
function  to  be  approximated,  the  structure  of  the  Galerkin  equations,  and  the  possibility 
to  use  fast  Fourier  transform  methods  when  evaluating  certain  integrals,  gives  an 
asymptotic  operation  count  for  determining  and  solving  the  Galerkin  equations  of  0( n 
log(n))  operations,  as  n  ♦  where  n  denotes  the  number  of  basis  functions. 

The  paper  is  organized  as  fallows.  In  Section  2,  we  present  the  families  of  near¬ 
orthonormal  functions  when  D  is  a  region,  whose  complement  is  bounded  and  simply 
connected.  Section  3  discusses  the  case  when  D  is  bounded  and  simply  connected.  In 
Section  4  the  Galerkin  scheme  is  presented  for  simply  connected  regions.  Section  5 
generalizes  the  scheme  to  multiply  connected  regions,  and  Section  6  contains  numerical 
examples.  An  appendix  contains  proofs  of  theorems  and  lemmas  of  Sections  2  and  3. 
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(2.5)  s(z)  2ir  ]  a(()|dc| 

by  Integration  along  3R  in  counter-clockwiaa  direction.  The  mapping  (2.5), 

s  ;  3(1  [0,21I[  has  an  inverse,  which  defines  a  parametric  representation  of  3(1.  If 

nothing  else  stated,  z(s),  0  <  s  <  2it,  refers  to  this  parametric  representation 

throughout  this  section.  The  domain  of  z(s)  is  extended  to  all  real  s  by  requiring 

z(s>  to  l>e  2z-perlodic.  It  is  convenient  to  also  consider  the  map  w 

z  e  3(1.  This  map  can  be  extended  to  (1^  as  follows.  Let  H(z,Wj)  denote  the  conjugate 

harmonic  function  to  G(z,Wj)  satisfying  H(Cg,  w^)  «  0.  Then 

♦  j(z)  exp(*K;(z,Wj)  -  iH(z,Wj))  maps  3(1  conformally  on  |w|  <  1,  is  continuous 

and  bljective  on  (1  ,  and  ♦.(w.)  =  0,  ♦.{?  )  =•  1.  The  sought  extension  of  w  « 
c  J  j  Jo 

is  given  by 

(2.6)  w  :-  (  n  z  e  (1  , 

j-1  ^ 


where  we  select  the  root  so  that  the  argument  of  w  varies  continuously  with  z  as 

‘  1/t 

traverses  3(1,  and  arg  (  11  ))  *  0. 

j-1  ^  ° 


Definition  2.1 

Let  m  be  a  non-negative  integer  and  6  a  real  constant  0  <  8  <  1.  C^'^ta.b) 

denotes  the  set  of  functions  that  are  m  times  continuously  differentiable  on  [a,b] 
with  the  mth  derivative  satisfying  a  Lipschitz  condition  with  Lipschltz  constant  6 
Let  z*(t),  0  <  t  <  L,  be  a  parametric  representation  of  3(1  w.r.t.  arc  length 
He  write  3(1  e  c“'®  if  z*(t)  e  c"'®tO,LI. 

The  following  lemma  is  needed  for  proving  some  of  the  subsequent  theorems. 

Lemma  2.1 

If  3(1  e  0  <  s  <  1,  n  integer  >  0,  then  z(s)  e  [0 ,2v  1  and 

o(z)  >  8  on  3(1  for  some  constant  8  >  0. 


z 


t. 


-3- 


The  next  theorem  Is  essential  for  the  description  of  the  properties  of  the  function 


families. 

Theorem  2.1 

Let  30  e  0  <  6  <  1,  n  integer  >  0.  Define 


(2.7a)  '■ 


n(2-s( — 5* - )) 

1-1 


n  (s-w,) 
k-1 


k/t 


k  -  1,2,3,. 


(2.7b) 


,,  ,*(2i-2).. 

n(z-i( i- — )) 

1-1 

0|j,2(*)  •“  4  ,  k  —  1,2,3,...  , 


n  (s-w  ) 
j-1  3 


k/l 


where  the  roots  are  selected  as  for  (2.1).  nien  for  some  constant  a,  0  <  a  <  2t, 
0^  ,(2(s))  •  e^''“(1+e"^''®)(1+0(k"""*))  as  k  ♦  uniformly  for  0  <  s 


(2.8) 


“  e^''*(1-e'^’'*)(1t0(k"""®))  as  k  uniformly  for  0  <  s 


For  *  8  0, 


(2.9) 


^,1(z)  “  «^'''*'*0<k*"“'”®)  as  k  ♦  *,  uniformly  for  z  in  a 

closed  subset  of  0, 

®k,2^**  -  e^*'**+0(k  "  ^  ®)  ag  k  ♦  **,  uniformly  for  z  in  a 

closed  subset  of  0. 


For  z  e  0  \30,  let  w  be  defined  by  (2.6).  Then 
c 


<  2«  , 

<  2z  . 


-4' 


(2.10) 


«  e^*'**(1+w  '')(1+0()c  "  ^  ®))  as  k  ♦  “/  uniformly  for  z  in  a 


closed  subset  of  0  \  3(1  . 

c 


®k,2**'  ”  e^*^*(1-w  *')(1+0(k  "  '  ®))  as  k  ♦  uniformly  for  z  in  a 


closed  subset  of  (1  \  3(1  . 

c 

If  z(s)  is  analytic,  then  the  convergence  in  (2.8)-(2.10)  is  geometric,  i.e.  there  is 


an  r,  0  <  r  <  1,  such  that  the  0(k  and  Otk  ”)  terms  can  be  replaced  by 


0(r’'). 


Remark  2.1 

The  constant  a  in  Theorem  2.1  depends  on  the  shape  of  (1  and  on  the  location  of 
the  Wj:s.  Prom  the  proof  of  the  theorem  it  follows  that  if  the  region,  the  Wj,  and 
z(0)  lie  symmetric  w.r.t.  the  x-axis  then  a  «  0.  ■ 

By  combining  the  Pk,i(z)  ^n6  CV,2***  various  ways  we  obtain  families  of 

functions  which  are  near-orthonormal  with  respect  to  the  inner  product 


2x 


(2.11) 


(f,g) 


z(8))g(z(s))  ds  , 


where  the  bar  denotes  complex  conjugation.  The  next  theorem  provides  examples  of  families 
of  type  (2.1).  The  theorem  follows  Immediately  from  Theorem  2.1. 

Theorem  2.2 

Define  the  families 

qo(z)  :=  1 

qy(z)  Qy^.,(z)  -  e  ,  k  “  1,2,3,...  , 


(2.12) 


and 


(2.13) 


qo(*)  1 

q^(z)  V2  (Ox,l(z)  -  k»  1,2,3,...  . 

For  members  of  both  of  families  (2.12)  and  (2.13),  we  have  under  the  same  smoothness 
recnilrements  on  z(s)  as  in  Theorem  2.1  that 


.1 

I 
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(2.14a)  q|^(z(s))  -  e  =  0(k  "  ®),  k  ♦  uniformly  for  0  <  s  <  2ii  , 

*’*11“  1  “fi 

(2.14b)  qk<*)  “  '^)  .  k  +  «"i  for  z  e  0,  uniformly  for  z  in  a  closed  subset 


of  £!  . 


(2.14c)  q)t(z)  “  e^''*w''(  1+0(k  "  '  ®)),  k  ♦  for  z  e  uniformly  for  z  in  a 


closed  subset  of  0  \  3(2  . 

c 

w  =  w(z)  is  defined  by  (2.6) 


If  z(s)  is  analytic,  then  the  0(k  ”  ® )  and  0(k  ^  ® )  terras  can  be  replaced  by 

0(r'')  terms,  where  r  is  a  constant  0  <  r  <  1. 

Example  2.1.  Let  (2  :=  {z:  |z|  >  1},  t  =  1,  and  w.  =  0.  Let  z(0)  =  1. 

Then  >()(zl  =  z,  i.e.  z(s)  =  e^®  and 


Qk,l(*)  !■= 


0k.2<^> 


„Cz-s(Il^)) 

- 

z  z 


n  (z-z(Ii|^)  , 

i=i  2k  ^k., 


=  1  -  z"’'  . 


By  symmetry  a  =  0,  and  from  (2.12) 


By  (2.13) 


q^(z)  i=  1  +  z-*^  -  1  =  z-’'.  k  >  1 


q^(z)  t=  V2  ( 1  +  z"’'  -  ( 1  -  z"*'))  =  z"*',  k  >  1 


I 


Hence,  in  this  case  the  families  of  functions  defined  by  (2.12)  and  by  (2.13)  agree. 


As  is  seen  ..n  the  proof  of  Theorem  2.1,  the  ordo-terms  in  Theorem  2,1  originate  from 
the  substitution  of  an  integral  by  its  trapezoidal  sum.  We  therefore  call  the  ordo-terms 
in  (2.8-10)  truncation  error.  In  the  members  of  family  (2,13)  the  truncation  errors 
cancel  in  part.  This  can  be  shown  by  studying  error  t>ounds  for  the  truncation  error  with 
the  Euler-MacLaurin  summation  formula,  see  Henrici  (1).  It  is  also  Illustrated  in  the 
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next  example.  The  error  is  measured  in  the  maximum  norm 

(2.15)  lf(z(8))<^  !“  sup  |f(z(s))|  . 

0<s<2z 

2 

Example  2.2.  Let  0  :•  {z«  x+  iy,  (^)  +  >  l}.  Choose  £  •  3,  and  let  “  0, 

wj  “  1.5,  and  Wj  «  -1.5.  Let  z(0)  =  2.  Then  by  synnetry  a  «  0. 

a)  defined  by  (2.12). 

-iXs 

k  Iq|^(z(s))  -  e 

1  6.8  •  10"’ 

2  2.3  •  10"’ 

3  7.7  •  10”^ 

4  2.6  •  10"’ 

5  2.7  •  10"’ 

10  2.2  •  10"^ 

20  5.2  •  10"^ 

30  3.5  •  10*® 

b)  defined  by  (2.13). 

-IVa 

k  lq^(z(s))  -  e  . 

1  3.5  •  10*'' 

2  1.4  •  10*' 

3  7.7  •  10*^ 

4  1.3  •  10*' 

5  1.4  •  10*' 

10  1.2  •  10*^ 

20  2.6  •  10*^ 

30  1.8  •  10*® 

The  cancellation  of  truncation  error  in  (2.13)  is  clearly  visibly.  ■ 

Henceforth,  ^)((z)  will  denote  functions  defined  by  (2.13).  We  next  introduce  a 
family  of  functions  of  type  (2.2).  Let 
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1, 


*  0 0  ^  $2  f  • » • 


(2.16) 


t 

n  (3 
j"i 


.k/t 


“»  k  *  1(1)t“1,  n  *  f 


n  (z-w  ) 

J-1  3 

where  the  root  is  selected  as  in  (2.13).  Then 

(2.17)  “  0/'«2,... 

is  an  example  of  a  family  of  type  (2.2).  Finally  we  introduce  a  family  of  type  (2.3)s 

[  V*)  1  . 


(2.18) 


^k'*> 


®k  l‘*’  ■  ®k  2*** 

k  M  1  7  t 

®k,l‘*’  *  ®k.2‘*’  '  '  . . 


We  conclude  this  section  stating  some  properties  of  the  function  families  Introduced.  Me 
begin  with  the  functions  (2.13). 

Theorem  2.3 

a)  If  for  all  k  >  k^,  lq|^(z(s))  -  e  then  ^ “Ijt ^ ^ ^ ***  linearly 

o 

independent,  and  the  numerators  of  the  OijSS,  k  >  k^,  have  all  zeros  in  0. 

b)  The  set  {q|^(z)}Q  is  linearly  dependent  only  if  for  aone  k 

^  k/t 

qv(z)  •  (z-w, )  vanishes  atz«Wj,j“1(1)t.  e 

J=1  J  ^ 

The  next  theorem  justifies  the  Galerkln  scheme  described  in  Section  4. 

Theorem  2.4 

Let  A(n)  be  the  linear  space  of  functions  analytic  in  0,  continuous  on  0,  and 

—  3  8 

with  piecewise  continuous  derivatives  on  0.  Assume  that  z(s)  8  C  '  [0,2ii], 

0  <  8  <  1.  For  f(z)  e  A(fl)  define 


2x 


J  f(z(s))e  ^*'*ds 


0.  1,  2, 


Then 


(  1  ®  norm  in  A(n). 

k-0  ^ 
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Remark  2.2 


We  only  have  to  establish  that  “  0  implies  f(z)  =0  on  Note  that 

0 

Parseval's  equality  might  not  hold;  f(z(s))  can  have  nonvanishing  Fourier  coefficients 
with  negative  index.  The  smoothness  requirements  of  f(z)  on  0  have  been  chosen  for 
convenience  and  can  be  relaxed.  a 

Theorem  2.5 

If  the  set  linearly  independent  and  if  r{s)  satisfies  the 

smoothness  requirements  of  Theorem  2.4,  then  the  function  set  ^ complete  in 

A(n)  equipped  with  norm  (2.15).  a 

We  next  turn  to  the  rational  functions  with  fixed  poles  defined  by  (2.17). 

The  q^^  are  near-orthonormal  in  the  sense  that  for  any  £  >  0  there  is  a  k^  such  that 

(2.19)  |(q^,  q.)l  <  £  Vk,j  s.t.  fk  -  j|  >  kj,. 

Moreover  ♦  0,  k  ♦  for  z  e  8,  uniformly  for  z  in  a  closed  subset  at  8. 

We  conclude  this  section  with  some  properties  of  the  rational  functions  ( z) 
defined  by  (2.18). 

Immediately  from  Theorem  2. t  we  obtain 
Theorem  2.6 

Let  z(s)  be  as  in  Theorem  2.1.  Then 

—4  Ve  -n“6 

q^(z(s))  =  e  +  0(k  )  #  k  uniformly  for  0  <  s  <  2ti  , 

q^(z)  =  0(k  )  ,  k  ♦  z  €  n,  uniformly  for  z  in  a  closed  subset  of  , 

qj^(z)  -  +  0(k  ^  ^  ^))#  k  ♦  *f  z  P  uniformly  for  z  in  a  closed 

subset  of  Q  \3fl,  w  =  w(z)  is  defined  by  (2*6)  • 
c 

If  z(s)  is  analytic#  then  there  is  a  constant  r#  0  <  r  <  1  such  that  the 
0(k  ^  and  0(k  ”  terms  can  be  replaced  by  0(r*')  terms.  ■ 
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Remark  2.3 


When  using  the  check  that  all  poles  of  the 

By  analogous  results  for  the  functions  there  is  a  kQ  depending  on  the  shape  and 

smoothness  of  30  such  that  for  k  >  kg,  all  poles  of  0^\30.  ■ 

The  only  depend  on  the  parametric  representation  of  30,  but  not  explicitly 

on  the  Wj.  As  the  next  theorem  shows,  we  can  choose  a  parametric  representation  without 
making  use  of  points  Wj . 

Theorem  2.7 

Let  z(s),  0  <  s  <  2x,  be  a  parametric  representation  of  30,  extend  z(s)  to 
a  2w-periodic  function  and  assume  that  z(s)  G  ®  (0 , 2it ]  ,  0  <  s  <  1,  n  integer  >  0, 

and  that  !z'(sll  >  0,  0  <  s  <  2ff.  q^^Cz)/  defined  by  substituting  the  present  parametric 
representation  z(s)  into  (2«18),  satisfies 

(2.19)  q^(z(s))  =  e“^^®(1  +  0<k"^“®)),  k  ♦  ®,  uniformly  for  0  <  s  <  211  , 

If  z(s)  is  analytic,  let  S  be  a  strip  (z  :  |lm{z>|  <  5)  on  which  z(8>  is 
analytic.  There  is  a  constant  r,  0  <  r  <  1#  depending  on  6,  such  that 


(2.20) 


qj^(z(s))  *  +  0(r*')),  k  + 


uniformly  for  0  <  s  <  2ir . 


3.  Bounded  simply  connected  regions 


We  now  consider  approximation  of  functions  which  are  analytic  on  bounded  simply 
connected  open  regions  and  smooth  on  !*  where  3S1  denotes  the  smooth 

boundary  of  !1^ .  Near-orthonormal  functions  are  derived  that  are  analogous  to  those  of 
Section  2.  The  special  case  when  the  near-orthonormal  family  consists  of  polynomials  has 
been  discussed  previously  in  [21.  The  proofs  of  the  theorems  are  similar  to  those  of 
Section  2,  and  moat  of  them  are  omitted. 

Let  To  define  the  function  families  select  points  Wj,  j  «  1(1)1,  in 

n  \3n.  Let  G(z,v^)  denote  the  Green's  function  for  the  Laplace  operator  for  with 

e  J  ® 

a  logarithmic  singularity  at  Wj,  i.e.,  &G(z,Wj)  »  0  for  t  6  {w^}  ,  G(z,Wj)  “  0 

for  z  e  3(1,  and  J  4^  (J!,w^)|d2|  -  1,  where  denotes  the  normal  derivative  into 

3(1  J 

(1  .  Then 
e 

1  ‘  3G 

(3.1)  a(z)  :•  J  I  (r,w^),  z  e  3(i  , 


Let  $  be  an  arbitrary  but  fixed  point  on  311,  and  define  for  z  B  3(1,  s(z)  by  (2.5) 
o 

with  0  defined  by  (3.1).  s(*)  has  an  Inverse  on  3(1,  and  we  denote  this  Inverse  by 

z(8)  1  [0,2x(  ♦  3(1.  We  define  z(s)  for  all  real  s  by  requiring  z(s)  to  be 

2x-perlodlc.  This  Is  the  parametric  representation  to  be  used  In  this  section.  Let 

H(z,Wj)  denote  the  harmonic  conjugate  to  G(z,Wj)  satisfying  HlS^^.w^)  -  0.  Then 

?j(z)  exp(-G(z,Wj)  -  1  H(t,Wj))  is  the  conformal  mapping  from  (1^  onto  |wl  >  1 

such  that  f  (w.)  -  ",  T  (C  )  -  1.  An  extension  of  w  »  to  (1  is  defined  by 

1  j  j  o  e 

(3.2)  w  :-  (  n  ?  (z))’''*,  z  e  (I  , 

J-1  ’ 


where  ve  select  the  root  so  that  the  argument  varies  continuously  with  z  z(s), 

*■  1/t 

0  <  s  <  2x,  and  arg  (  H  T  (E  ))  "0.  The  next  theorem  parallels  Theorem  2.1. 

j-1  3  ° 
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Theorem  3. 1 


Let  z(3)  e  0  <  P  <  1,  n  integer  >  0.  Define 


(3.3a) 


== 


i«i 

C*'  f  (3-.)’'^* 
n  3 
j“i 


h  -  1,2,3,. 


(3.3b) 


0k,2(z)  :■ 


n,.-z(ll|f^), 
j  =  1 


k  -  1,2,3,. 


The  branches  of  the  roots  are  selected  so  that  the  arguments  vary  continuously  as  z 
traverses  3(1.  c  denotes  the  capacity  of  3(1  relative  to  point  charges  at  z  •  Wj, 
j  =  1(1)1,  i.e.  c  is  part  of  the  solution  {c,u)  of  the  system  of  integral  equations 


(3.4) 


In  —  +  J  ln|z-C|u(5)|<15l  =  -r  J.  ln|z  -  w  |,  z  e  3(1  , 
3S1  *  j»=1  ’ 


j  u(C)|d«;l  -  1. 
3(1 


Then  for  some  constant  a,  0  <  a  <  2t, 

C*ir  i(z(s))  =  e^''**(e^*'*  +  1)(1  0(l5  "  ®)),  It  ♦  •,  uniformly  for  0  <  s  <  2* 

(3.  5)  •  '  IV  ik  fl 

*  e^''“(e^  -  1)(1  +  0(k"""‘’)),  k  ♦  uniformly  for  0  <  s  <  2it 

For  z  e  (Ij^ , 

(?V  i(*)  =  e^*"*  +  0(k  "  '  ®),  k  ♦  uniformly  for  z  in  a  closed  subset  of  (1, 
(3,6)  '  ilc  -  -1-6 

Q|j  2^*1  “  “S  **  +  0(k  "  **),  ){-*•",  uniformly  for  z  in  a  closed  subset  of  (1^ 

For  z  e  (1  \3(l,  let  w  be  defined  by  (3.2).  "nien 
e 
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1 


Ojj  ^(z)  "  +  1)(1  +  0<k  ”  '  ®)),  k  ♦  unifomiXy  for  z  In  a 

closed  subset  of  Sn  . 

®k  2**'  “  “  1){1  +  0(k  **  ^  ^))i  k  ♦  *,  uniformly  for  z  in  a 

closed  subset  of  S)  N  3fl  . 

e 

If  z(s)  Is  analytic,  then  the  convergence  in  (3.S)-(3>7)  Is  geometric,  i.e.  there  Is 
constant  r,  0  <  r  <  1,  such  that  the  0(k  ’^)  and  0(k  )  terns  can  be  replaced 

by  0(r’')  terms.  ■ 

Similarly  as  in  Section  2,  vre  can  now  combine  the  knd  various 

ways  to  obtain  families  of  near-orthonormal  functions,  with  respect  to  an  inner  product 
(2.11)  with  z(a)  defined  as  the  Inverse  of  s(z). 

Theorem  3.2 

Define  the  family 

qo(*)  1 

q|j(z)  :•  *2k,2***  "  ^»2,3,...  . 

Under  the  smoothness  conditions  of  Theorem  3.1, 

(3.9a)  q^(z(s))  -  e^''*e^''®  +  0(k'"*®),  k  *  ",  uniformly  for  0  <  s  <  2ii  , 

1 •fi 

(3.9b)  <Ik^*'  “  ) ,  k  ♦  z  e  0^,  uniformly  for  z  in  a  closed  subset  of 

{3.9c)  +  0(k”"*'”®)),  k  »  z  e  n^\3n,  uniformly  for  z  in  a 

closed  subset  of  0  \  3(1.  w  •  w(z)  is  defined  by  (3.2)  . 
e 

If  z(8)  is  analytic  then  the  0(k”"”®)  and  0(k  "  ’  ® )  terms  can  be  replaced  by 
0(r'')  terms,  where  r  is  a  constant  0  <  r  <  1.  ■ 

Theorem  3.3 

a)  If  for  all  k  >  k  ,  lqj^(z(8))  -  <  1,  then  ^'’x’^'^k-k  linearly 

o 

independent,  and  the  numerators  of  the  qj^  have  all  their  zeros  in  (1^,  X  *  X^ . 

b)  The  set  {q  (z))"  Is  linearly  dependent  only  if  for  some  k 

f  k  0 

-  X/i 

qv(t)  *  (z-w  )  vanishes  at  z  »  w, ,  j  “  1(1)1. 

j-1  3 
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Theorem  3.4 


Let  f(z)  e  A(n^),  where  is  defined  analogously  to  A(n},  c.f.  Theorem 

3  d 

2.4.  Assume  that  z(s)  6  C  '  t0,2itl,  0  <  B  <  1.  Define 


1 

a^  J-  3^  i  f(z{8))el'‘®ds.  k-0,1,2 . 


Then  (  la  |^)‘^  Is  a  norm  In  A(n  ).  a 

k-l  ^ 

Theorem  3.5 

If  the  set  {q|^(z)}Q  is  linearly  Independent,  and  If  z(s)  satisfies  the  same 
smoothness  conditions  as  in  Theorem  3.4,  then  the  {qj^(z))p  are  complete  in  A(flj  ) 
equipped  with  norm  (2.15),  s 

We  next  turn  to  rational  functions  with  fixed  poles  corresponding  to  functions 
(2.17).  With  the  Wj.-s  of  the  present  section,  we  define  functions  hj(z)  as  in  (2.16), 
where  we  choose  the  root  so  that  the 

(3.11)  '•){’*’  :»  qj^izlh^Cz),  )c-  0,1,2,... 

are  rational  functions.  The  neat-orthonorraal  in  the  sense  of  (2.19).  Also 

q|^(z)  ♦  0,  l<  •*  “  for  z  e  uniformly  for  z  in  a  closed  subset  of  £5^^. 

Finally,  we  define  functions 

[  q(,(*)  1  . 


(3.12) 


q^(*) 


®K.1<*>-°k,2‘*’  ' 


1,2,3,. 


The  functions  (3.12)  are  obtained  by  inverting  the  functions  (2.18),  and  this  determines 
their  properties. 


-14- 


.  A  Galerkln  scheme  for  simply  connected  regions. 

We  describe  a  Galerkln  method  for  numerical  approximation  of  analytic  functions  on 
regions  il,  whose  complement  In  bounded  and  simply  connected,  as  In  Section  2.  The 
structure  of  the  resulting  matrix  makes  fast  numerical  solution  of  the  Galerkln  equations 
possible.  We  assume  that  the  parametric  representation  z(s)  Is  known.  A  method  for 
computing  z(s)  for  given  Wjts  Is  described  In  [3],  [4].  The  presentation  of  the 
approximation  method  for  analytic  functions  Is  followed  by  a  description  of  the  method 
applied  to  the  approximation  of  harmonic  functions.  We  begin  with  a  scheme  for  the 
approximation  of  functions  f<z)  e  A(fl)  by  a  linear  combination  of  qj(z):s  defined  by 
(2.13).  In  the  Galerkln  scheme  use  the  functions  1  ,e^*,  •  •  • /e^***”’ as  test 
functions  and  an  Inner-product  obtained  by  replacing  the  Integral  by  a  trapezoidal  sum 
with  N  nodes.  The  Galerkln  equations  become 
(4.1)  ^  ' 

where  A  «  ^^Jk^kij-O'  —  ” 

.  <  ,  ,  A  A  2* 

(4.2a)  V'-TTn  i  qj(*(mis))e  .  &•  :- 

•'  ra“0  ■' 

N-1  _  . 

(4.2b)  b  I-  I  f(z(mAs))e  *  . 

maO 

We  consider  the  case  when  f(z(s))  Is  such  that  a  large  system  (4.1)  has  to  be 
solved  In  order  to  obtain  an  accurate  approximation  and  show  how  (4.1)  can  be  set  up  and 
solved  In  0(N  log(N))  multiplications  or  divisions.  Let  c  >  0  be  arbitrary  and 
partition  the  N  *  N  matrix  A  as 

*11  *12 
*21  *22 

where  A^i  Is  an  x  Ni  matrix  with  choosen  so  that 

maxdl-Ajjl,  IA,2*>  ’*21*^  *  ^  ' 

where  I  I  denotes  the  maximum  norm.  For  smooth  contours  we  can  select  Independent 
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of  N.  A 


^2  and  A2^<  we  replace  by  zero-matrices,  and  A22  we  replace  by  the  identity 
matrix.  This  modified  system  can  be  solved  in  O(H^)  operations.  The  elements  of 
and  of  _b,  we  determine  by  a  fast  Fourier  transform  method.  This  requires  0(M  logCN)) 
multiplications. 

The  matrix  A^^  can  be  ill-conditioned.  In  those  cases  the  modification  of  the 
system  can  be  regarded  as  a  regularization:  it  prevents  the  coefficients  of  the  q^, 
k  >  to  grow  with  N. 

The  operation  count  for  solving  the  subsystem  with  matrix  A^ ^  can  be  reduced  by 
using  an  iterative  solution  method.  Iterative  methods  may  be  advantageous  since  Af|  has 
the  same  structure  as  A.  For  a  discussion  on  suitable  Iterative  schemes,  see  [31. 

The  above  discussion  Is  valid  also  when  we  seek  an  approximant  in  the  fora  of  a 
linear  combination  of  function  q|^(z).  When  functions  are  used,  we  write 

A22  ”  B22  ♦  C22»  where  'Cjjlp  <  £  and  B22  *-8  a  banded  matrix  of  band  width 

independent  of  N.  Ihe  operation  count  therefore  also  in  this  case  is  0(N  log(N)). 

Next  we  consider  the  approximation  of  harmonic  functions  g  in  0  by  a  linear 
combination 


N  N-1 

C)(z)  :•  Cq  +  I  Cj),  Re(q^(z))  +  Ira(q^(z))  . 

k-t  k-1 


We  require  for 


w 


1. 


(4.3) 


2N-1 

I  (0(r(Sj))  -  g( z( s j ) ) )cos( ms j )  *0,  m  •  0(1)N 
1*0 

2N-1 

I  (0(z(s.))  -  g(z(s.) ) )sin(ms^)  -0,  m  -  1 ( 1 )N  -  1  . 
J-0 


These  Galerkin  equations  can  be  ordered  to  have  the  same  structure  as  (4.1)  and  the 


operation  count  is  0(N  log(N))  for  setting  up  and  solving  (4.3). 


5.  A  Galerkln  scheme  for  multiply  connected  reqlona. 


M 

Let  SI  have  simply  connected  complements  with  smooth  boundary  curves 

3Slj^.  A  function  f(z)  analytic  In  n,  can  be  partitioned  into  M  parts 


f^(t) 


1 

211 1 


liLL 

5-z 


dC, 


z  e 


complement  of 


k 


1(t)M  , 


with 


M 

f(z)  »  \  *k^**  °0'  *  ®  *  constant  . 

k=*1 

Each  of  the  can  be  approximated  by  the  functions  of  sections  2  or  3.  In  the 

Galerkln  scheme  for  approximating  f(z),  we  use  basis  functions  for  the  complement  of 
Slj^,  k  •  1(1  )M.  Since  the  functions  g^lz)  for  the  complement  of  SJ^^  tend  to  0  as 
i  on  each  of  the  boundaries  3S1  ,  m  +  k,  the  Galerkln  equations  can  be  ordered  so 

that  the  same  structure  as  for  system  (4.1)  is  obtained.  Approximation  of  harmonic 
functions  is  carried  out  analogously,  but  we  note  the  need  for  functions  with  logarithmic 
singularities,  see  Walsh  [6]  for  details. 
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6 .  Numerical  examples* 

The  examples  illustrate  some  of  the  properties  of  the  functions  and  q^^vz). 

All  computations  were  done  on  a  VAX  11/760  in  double  precision,  i*e.  with  12  significant 
digits. 


Example  6.1. 

Let  w . 

3 

and  z(s)  be 

the  same  as  in  example  2.2,  see  Figure  6.1 

Q 

o\ 

Figure  6.1 

Let 

(6.1) 

2_ 

^kk  2. 

2it 

J  q^(z(B))e"^*'®ds  , 

0 

(6.2) 

2n 

^k  ••=  ^  i 

0 

k 

^kk 

1 

1.33 

2.70  •  10"^ 

2 

1.11 

7.51  •  10'“ 

3 

1 .04 

1.50  •  10"^ 

4 

1.02 

5.53  •  10"^ 

5 

1.03 

5.53  •  10"^ 

10 

1  +  1.3  • 

10"^  2.24  •  19"^ 

20 

1  +  7.6  • 

10"®  3.29  •  10*® 

30 

1  +  5.8  • 

10~’°  2.94  •  10"’° 
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The  table  indlcatea  that  all  off-dlaqonal  elemante  In  the  Galerkln  equations  are 
small,  and  that  Iterative  solution  can  be  employed  for  all  slses  of  the  linear  system 
(4.1). 


Example  6.2.  Let  0  be  the  complement  of  a  region  formed  by  a  rectangle  with  sides  of 
length  2  and  1  and  2  disc  halves,  see  Figure  6.2.  Let  i  -  3  and  >  0,  W2  ~  1.5, 

W3  -  -1.5,  z(0)  >2.  By  symmetry  a  -  0  In  Theorem  2.1. 

n 

•►Z(O) 


Figure  6.2 


Let  and  dj^ 


be  defined  as  In  (6.1)-(6.2). 


k 

lq^(z(s)  )-e"'‘’‘*l. 

*kk 

Ok 

1 

3.50 

lO"'* 

1.27 

4.74 

•  10 

2 

1.39 

10"'’ 

1.01 

1.20 

•  10 

3 

2.44 

10"’ 

1.01 

3.33 

•  10 

4 

3.83 

10“'' 

1.07 

5.59 

•  10 

5 

9.29 

10“2 

9.92  •  10“' 

3.41 

•  10 

10 

4.64 

10“2 

1 

+  4.9  •  10"^ 

9.97 

•  10 

20 

2.44 

0 

1 

1 

♦  8  •  10"^ 

0 

0 

•  10 

30 

1.60 

10“2 

1 

-  2  •  10“^ 

3.07 

•  10 

The  convergence  of  qy(z(s))  to  Is  alow,  however,  the  Galerkln  equations  are 

strongly  diagonally  dominant  already  for  small  systems  of  equations.  ■ 
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1 


Example  6.3.  Let  St  be  the  exterior  of  the  curve 

iSl  :»  {z  •  x+iy,  x  “  1.4co8(t)  +  2.8coB(2t)  -  2.8, 

y  -  2.45sin(t-0.2)  +  1.66in(t)  +  0.98Bin(2t)  -  0.49Bin(4t)  -  0.74, 

0  <  t  <  2x}  , 

Bee  Figure  6.3.  t  »  5  and  the  Wjta  are  (-4,4),  (-1.8, -2. 9),  (0,-1),  (-1,1),  (-2.4, 
2.2)  and  (-4.5,3). 


For  this  example  a  ^  0  In  Theorem  2.1.  Instead  of  determining  a  and  showing  a 
table  for  z(  s)  )-e^*"*e  we  have  determined,  by  the  argument  principle,  the 

number  of  zeros  of  the  numerator  of  St. 
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k 


*  of  zeros  of  numerator  In  SI 


!*kkl 


1 

5.49 

10"’ 

6.66 

10“* 

0 

2 

1 

- 

5.09 

10"2 

4.04 

10"’ 

0 

3 

1 

+ 

2.08 

10“’ 

2.15 

10"2 

0 

4 

1 

' 

2.11 

io~2 

1.21 

10"’ 

0 

5 

1 

■V 

6.28 

10~2 

1,12 

10"’ 

0 

10 

1 

- 

1.2 

10“^ 

1.S0 

O 

1 

K> 

0 

20 

1 

- 

3.4 

O 

1 

1.15 

10"^ 

0 

30 

1 

5.4 

10“^ 

5.74 

lO"* 

0 

■ 

Example  6.4.  We  consider  approximation  on  a  bounded  region, 

0^  s«  {z  »  X  t  iy,  bet  the  Wj  »  j  -  3  +  31,  see  Figure  6.4.  This 

allocation  of  the  v^:a  is  reasonable  If  the  function  to  be  approximated  is  known  to  have 
Its  closest  singularities  to  3S1  In  the  vicinity  of  the  Wjjs.  Again  a  In  Theorem  2,1 
Is  nonzero.  Kow 

2v 

(6.3)  *kk  *“  27  /  qk(z(s))e^’'*d8 

0 

1 

(6.4)  ^k  ?7  I  *  l»kk'^ 

+  +  4-  +  + 

w. 

3 


Figure  6.4 


k 


l^kkl  ''k  »  of  zeros  of  numerator  of  i*' 


1 

1 

-  4.45 

10"^ 

3.50 

10“^ 

1 

2 

1 

-  3.22 

10"^ 

1.12 

10'** 

2 

-3 

-2 

3 

1 

-  5.32 

10  ■* 

4.08 

10  ^ 

3 

4 

1 

-  2.78 

10'^ 

1.50 

10-2 

4 

5 

1 

-  3.85 

10“^ 

7.84 

10-2 

5 

10 

1 

-  5.63 

10-« 

1.13 

10“^ 

10 

20 

1 

-  2 

10'^^ 

2.35 

10-« 

20 

30 

1 

-  2 

10-1^ 

5 

10-’2 

30 

The  remaininq  examples  treat  the  functions  g^^Ct). 

Example  6.5.  Let  fi,  i,  wj  and  z(s)  be  the  same  as  in  example  6.1.  Tlie  aj^j^  and 
djj  are  defined  by  (6.1)-(6.2)  with  q^(z(s))  replaced  by  q^(z(s)). 


k  #  of  zeros  in  fl  of  ^*kkl  “^k  lq^(z(s))-e 

numerator  denominator 


1 

1 

0 

1.34 

2.22 

lO"’’ 

1  -  3 

10"^ 

2 

2 

0 

1.11 

1 .69 

io"^ 

2.59 

10“*’ 

3 

3 

0 

1 

+ 

4.0 

• 

10-2 

1.47 

10“^ 

7.69 

10"^ 

4 

4 

0 

1 

+ 

6.6 

• 

10-2 

1.13 

10*^ 

1.69 

10“’ 

5 

5 

0 

1 

+ 

6.3 

• 

10-2 

2.  10 

TO'^ 

2.  70 

io“’ 

10 

10 

0 

1 

+ 

2.6 

• 

10-^ 

1.32 

10*^ 

2.40 

10"^ 

20 

20 

0 

1 

•I- 

6.6 

• 

io-« 

6.58 

10'» 

5.09 

lO'"* 

30 

30 

0 

1 

7.8 

• 

lo-’o 

5.87 

10-^'^ 

3.53 

in'^ 

22- 


^cr 


Example  6.6.  Let  and  z(a)  be  the  same  as  in  example  6.2. 


k 

#  of  zeros  in  Q  of 
numerator  denominator 

u 

kk 

■Jk 

lq|^(z<s))-e"^''® 

1 

1 

0 

1.31 

1,90 

•  10"’ 

1 

u> 

o 

t 

2 

2 

0 

1  - 

1.4 

• 

10“2 

2.01 

•  10'^ 

2.32  •  10"’ 

3 

3 

0 

1  - 

8.9 

• 

10"^ 

5.71 

1 

o 

• 

3.75  •  10"’ 

4 

4 

0 

1  - 

3.4 

• 

10"2 

9.89 

•  10“^ 

5,42  •  10"’ 

5 

5 

0 

1  - 

7.8 

• 

10"2 

7.01 

•  10"^ 

1.63  •  10"’ 

10 

10 

0 

1  * 

2.8 

• 

10"^ 

2.04 

1 

O 

• 

8.24  •  10"^ 

20 

20 

0 

1  + 

2.3 

• 

10*^ 

1.92 

•  10"* 

4.05  •  10'^ 

30 

30 

0 

1  - 

3.0 

• 

10~* 

5.90 

1 

O 

2.66  •  10"^ 

Example 

6,7.  Let 

fl,  1,  and 

z(fl) 

be 

the  same  as 

in  example  6.3. 

k  #  of  roots  in  fl  of 

numerator  denominator 


1 

1 

0 

6.06 

10"’ 

4.15 

• 

10"’ 

2 

2 

1 

7.99 

10"’ 

5.35 

3 

3 

0 

1 

+ 

8.90 

1 

o 

1.97 

• 

10"’ 

4 

4 

0 

1 

- 

1,43 

10"2 

3.76 

• 

10"’ 

5 

5 

0 

1 

2.80 

10"^ 

1.87 

• 

10"’ 

10 

10 

0 

1 

+ 

4.45 

10"^ 

3.21 

• 

10"^ 

20 

20 

0 

1 

+ 

2.81 

10"* 

1.84 

• 

10"^ 

30 

30 

0 

1 

• 

3.07 

10"5 

1.05 

• 

10"* 

^(z)  has  to  be  replaced,  since  it  is  singular  in  (2.  a 
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Example  6.8.  Let  12^,  t,  and  z(s)  be  the  sane  as  in  example  6.4. 


k  #  of  zeros  in  of  l*kkl  ‘^k 

numerator  denominator 


1 

1 

0 

8.79 

10"’ 

4.52 

• 

10-1 

2 

2 

0 

1 

+ 

4.45 

10"2 

3.33 

• 

10-1 

3 

3 

0 

1 

♦ 

1.74 

10"2 

8.39 

• 

10-2 

4 

4 

0 

1 

- 

1.00 

10"2 

2.92 

• 

10-2 

5 

5 

0 

1 

+ 

6.66 

10“^ 

1.59 

• 

10-2 

10 

10 

0 

1 

- 

5.61 

10“5 

2.26 

• 

O 

1 

20 

20 

0 

1 

+ 

1.59 

10"^ 

4.70 

• 

OD 

1 

O 

30 

30 

0 

1 

+ 

7.45 

10-12 

1.07 

• 

10-11 

The  concluding  examples  illustrate  that  z(s)  can  be  chosen  quite  arbitrarily  when 
approximating  by  functions  qj^(z). 
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Example  6.9.  Let  !)  be  the  region  of  example  6.2.  and  let  z(s)  be  the  arc  length  of 
3n  with  z(0)  >  2. 


k  #  of  ceroe  In  n  of 
numerator  denominator 


1  1  0 

2  2  0 

3  3  0 

4  4  0 

5  5  0 

10  10  0 

20  20  0 

30  30  0 


“kk  ^k 

1.36  1.48  • 
1  -  4.24  •  10"2  4.71  • 
1  -  9.23  •  10"3  1.11  • 
1  +  1.82  •  10"2  3,18  • 
1  -  3.51  •  10“^  4.98  • 
1  +  1.32  •  10"^  1.75  • 
1  +  5.26  •  10“^  2.50  • 
1  +  2.79  •  10“^  7.95  • 


10-’ 

1  -  4 

• 

10-4 

10~2 

3.03 

• 

10-'' 

10'2 

1.81 

• 

10-^ 

10"^ 

1.21 

• 

10-' 

10"^ 

3.36 

• 

10"2 

10-4 

3.65 

• 

10~2 

10-5 

2.04 

• 

10-2 

10-5 

1.36 

• 

10*2 

Example  6.10.  Let  id  and  z(8)  be  the  same  as  in  example  6.9,  but  compute  now 
fiuctlons  for  the  approximation  on  11^  «  the  Interior  of  3!i, 

k  #  of  zeros  in  of  l*kkl  "^k  z(  s) 

numerator  denominator 


1 

1 

0 

7.68 

10-’ 

3.93 

10*2 

5.00 

10*’ 

2 

2 

0 

1 

- 

1.56 

10*2 

4.93 

10-2 

3.12 

10"'' 

3 

3 

0 

1 

+ 

1.56 

10-2 

1.10 

10-2 

1.54 

10-' 

4 

4 

0 

t 

- 

1.85 

10*2 

2.76 

10*2 

1.11 

10-' 

5 

5 

0 

1 

+ 

3.62 

10-4 

5.00 

10-4 

3.48 

10-2 

10 

10 

0 

1 

- 

1.16 

10*2 

1.70 

1 

o 

3.52 

10-2 

20 

20 

0 

1 

- 

5.02 

10-4 

2.44 

10-5 

1.99 

10-2 

30 

30 

0 

1 

- 

3.73 

10-4 

7.83 

10*5 

1.35 

10*2 
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The  appendix  contains  proofs  of  moat  of  the  theorems  and  lemmas  of  sections  2 


Proof  of  Lemma  2 ■ i . 

The  correspondence  between  the  smoothness  of  the  ^^(z)  and  the  smoothness  of 
has  been  investigated  by  Warschawski.  See  Suetin  (SI  for  statements  of  results  and 
references,  n  >  0  implies  that  there  are  constants  5^  >  ^  such  that 

^2  ^  3  ”  KDl.  Since  s(z)  and  z(e)  are  as  smooth  as  the 

the  first  part  of  the  lemma  follows.  Prom  |♦•(z)l  «  2ii  {z,Wj).  z  S  3(1,  j  »  it 

follows  that  a(z)  >  S  on  30  for  some  constant  5  >  0.  s 

Proof  of  Theorem  2.t. 

Conformal  mappings  :  0^  *  (w|  <  1  with  "  **  satisfy  for  z  e  0^  and  for 

some  real  constant  ,  0  <  <  2it  , 


)  »  J  ln(— ^  (C.w  )|dc|  +  io., 


j  30 


^  are  the  same  functions  as  in  (2.4),  The  constants  depend  on  the  location 

of  Cg/  the  point  on  30  satisfying  ♦jCCg)  •  '»  j  “  Integrate  along  30 

from  Cg  and  choose  a  branch  of  the  logarithm  which  changes  continuously  with  the 

variable  of  integration.  Substituting  the  map  w  «  ♦(z)  ;«  z  into  (A.1)  yields 

(A. 2)  0  ^  i  |wl  <  1  » 

2”  |w|>1  w  -  i 

Average  the  1  equations  (A.l)  and  let  a  ^  •  By  (2.4)  and  (2.5)  we  have  for 

j-1  ^ 


ln(  n  ♦  (z)) 
j-1  J 


ln(  n  (Z  -  w  )  '  ) 
J-1  J 


Substitute  (2.6)  into  (A. 3),  and  subtract  (A. 2)  from  (A. 3).  We  obtain,  with  w  w(z). 


(A. 4) 


ln(w<z))  -  ln(  n  (z  -  w  )  j  ' 

j.1  ’  0  t  -  z(t) 


where  we  choose  a  continuous  branch  of  the  logarithm,  the  functions  And 

we  obtain  by  discretizing  the  integral  in  (A. 4)  by  the  trapezoidal  rule  and 
exponentiating.  We  Introduce  the  nodes 

tj  s-  (j  -  '1.  i  ‘  ’(’>N  , 


(A. 5) 


tj  ^  < j  -  j).  J  =  1(1)N  • 


The  discretization  error  depends  on  the  smoothness  of  the  integrand.  We  write 


(  A. 6) 


it 


z)e^''^,  z  e  tl 


Suetin  (5)  has  studied  the  derivative  w.r.t.  t  of  (A. 6)  for  £  •  1  in  (2.6).  Froei 

Lemma  2.1  and  Suetin  [S],  Lemma  1,  it  follows  that  if  £  •  1, 

(A. 7)  "  0(  |h  I  ,  |k|  ♦  uniformly  for  z  e  . 

Since  the  smoothness  of  z(t)  does  not  change  if  we  choose  a  finite  £  >  1,  (A. 7)  holds 

for  i  f  1  also.  When  z  e  fl  \3{1,  the  left  hand  side  of  (A. 6)  can  be  differentiated 

c 

once  more.  Therefore 

(A. 8)  ^  I*'!  *  •»  t  ®  0^\3B,  uniformly  for  z 

in  closed  subsets  of  (I  \  3(1  . 

c 

When  3(1  is  analytic,  we  have  for  some  constant  r,  0  <  r  <  1, 

(A. 9)  “)ct*)  ’  0(rl’'l),  |k|  ♦  *,  uniformly  for  z  6  0^  . 

We  next  utilize  the  connection  between  discrete  and  continuous  Fourier  coefficients 


(A. 10) 


-  I 

"l-l 


Ini 


w(z)  -  e 
z  -  z(t^) 


.(N) 


(z) 
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(A. 11) 


■o  ***'  *  i  *1cn(** 

)t— « 


We  obtain 

(A.  13a)  aQ**^(r)  -  aQ(i)  •  0(M  "  '  ®),  »)♦••,  uniformly  for  z  e  0^  , 

(A.  13b)  a^^’Nz)  -  a(,(z)  -  N  ♦  »,  z  e  n^\30,  uniformly  for  z 

in  closed  subsets  of  S)  \  30  . 

c 

Substituting  (A*10)f  (A*11)  and  (A*13a)  into  (A«4)f  multiplying  by  N  and  exponentiating 
yields 


N  it. 

J  n  (w(z)  -  e  ^) 

*(z)*»  n  (z  .  «  1-”/*  -  ^ - 

1-1  J  ” 

■'  n  (z  -  z(tj ) 

j-1  3 


e^““(1  t  0(N'''-®))  , 


N  ♦  »,  uniformly  for  z  6  30.  Now 


(A. 14) 


n  (w(z)  -  e  ^)  “  w(z)*^  •  1 
j-1 


n  (z  -  z(tj) 


0n,2<*>  “  I 


-  «^"“(1  -  w(z)-’')(1  ♦  0(N-"-®)).  N 


n  (z  -  w.) 

J-1  ’ 


This  shows  part  of  (2.8),  since  by  definition  w(z(s))  «  e^*. 


The  other  part  of  (2 .8)  we  obtain  by  uaing  the  nodes 

discretization  of  the  integral  in  (A.4).  When  z  8  VdQ,  we 

c 

and  (A. 13b)  into  (A»4)>  This  shows  (2.10).  To  show  (2.9),  we 

(A. 15)  ln(z  -  Wj)  -  j  ln(z  -  C)  (C»w  )|dc|,  j 

^  3fl  3n  j 


of  (A. 5)  in  the 
substitute  (A. 10),  (A. 11) 
note  that  for  zed, 

-  1(1)t  , 


and 


i  211 

(A. 16)  -  I  ln(z  -  wM  .  J  ln(z  -  C)<J(C)|di:|  -  -r-  J  ln( z  -  z( t )  )dt  , 

*  j-1  ^351  ^  0 


where  we  choose  continuous  branches  of  the  logarithms.  Now 


1  "  ikt 

—  J  ln(z  -  z(t))dt  *  ^  b|^(z)e  ,  zed, 

0  )c“-«» 


where  b)^(z)  =  0(  |  Ic  |  )  ,  |lc|  ♦  zed,  uniformly  for  z  in  closed  subsets  of  d. 

The  rightmost  integral  in  (A. 16)  we  discretize  using  the  trapezoidal  nodes  (A. 5)  and 
express  the  truncation  error  with  the  Fourier  coefficients  by(z),  similarly  as  in  the 
proof  of  (2.8).  This  gives  (2.9).  a 

Proof  of  Theorem  2.3. 

Part  a):  If  lq|^(z(s))-e  <  1,  then  the  winding  number  of  as  z  traverses 

3d  counter-cloc)5wise  is  -  )«.  This  shows  the  linear  independence.  Since  the  denominator 
of  q(t(z)  has  winding  number  -  k  aB  z  traverses  3d  counter-cloc)5wise  for  every 
k  >  0,  it  follows  that  the  winding  number  of  qj^{z),  k  >  k^,  is  that  of  the 
denominator,  and  therefore  the  numerator  cannot  have  any  zero  in  d  U  3d . 


Part  b) 


:  It  suffices  to  consider  linear  combinations  of  the  form  ^  a  ,  .  q  ,  .(z), 

t  ,  .  k=0 


where  j  is  fixed.  Let  Pj^^^  Jz)  •  H  (z-w^)''  Then 

m  n-1 


I  a. 


k=0 


ki+j^’ki+j 


(z)  SO  on  3d  is  equivalent  to 
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*  iR  k 

°"  *"•  "*  “y  assunia  that 

*m*+j*’mi+j'*'  ^  Then  has  zeros  at  z  “  Wj ,  j  »  1  ( 1  )l .  • 

In  the  following  proofs,  except  that  of  Theorem  3.1,  we  let  a  >  0,  since  this 
simplifies  the  formulas  somewhat  and  does  not  change  the  method  of  proof. 

Proof  of  Theorem  2.4 

We  have  to  show  that  there  Is  no  function  f(z)  in  A(0)  whose  Fourier  coefficients 
all  have  strictly  negative  Indices.  We  assume  the  contrary,  l.e.,  let  f(z)  e  A((i)  and 
assume  that  f(z(s))  has  the  representation 


(A.  17)  f(z(s))  •  I  a.i^e^*'®,  with  a.j^  ^  0  for  some  k  >  1  . 

k=-1 

For  each  Integer  N  >  1,  S|^(f(z))  e  Affl),  tdiere  s^  is  a  scaling  factor  choosen  so  that 
the  coefficients  al]J'  defined  by 


Sj,(f(z(s)))** 


I  al]J>ei’'» 

k-N 


satisfy  I  Define 

k»N 


c^(z) 


8„(f(Z))"- 

N 


I 

k”K 


•’k'*’ 


z  e  R  u  sn  . 


Then  both  Re(e^(z)),  Imltj^lz))  are  harmonic  In  fl,  piecewise  continuous  on  R,  and 
therefore  achieve  their  maximum  on  the  boundary.  By  selecting  N  sufficiently  large,  ve 
can  make  small  on  3R.  This  Is  seen  by  writing  t^Cz)  in  the  form 


e  (z(8) ) 
w 


I 

k-N 


(N)  Iks 
*-k  ® 


-  I 


(N) 


k-N 


)) 


k-N 


a'”’o(k"") 

-k 


By  the  smoothness  requirements  on  z(s),  we  have  n  >  1.  By  the  Schwarz  inequality 


V  k  >  1.  Hence 


|€  (Z(8))|^  <  I  i  0(k‘^")  -  I  0(k^")  . 

k-N  k-M  k-N 


There  la  a  constant  c  independent  of  k  such  that  0(k  <  c  k 


|e  <z(s)) 

N 


Select  N  such  that  •t„t*(a))l  <  - 

•*  2/1 

Bessel's  Inequality 


«• 

I  "2n. 

J  X  dx  »  c 

N 


2n-1  ■ 


implies  there  is  a 
is  a  z  e  n  where 


1 


I 

k-N 


l“.k  ' 


2x 

J 


0 


|s^,f^(z(s)  )  l^ds 


z  e  Jn  where  (sj,f**(z)|  >  1.  By  continuity  of  f**(z)  in  £1, 

|a„f"(z)|  >  ^6  •  Let  the  disk  | z-z I  <  S  lie  entirely  in  ft. 
N 


there 

Then 


(A. 18) 


1 

2xi 


f  (z) 


dz 


^  ^  i  al”’  q^(*)  )(*-*)■■' de 

A  k-N 

lz-z|-i 


1 

2xi 


I 

I  z-z 1-4 


dte 


N  ^ 

The  first  integral  on  the  left  hand  side  equals  s  f  (z),  the  second  vanishes.  The 

N 

integral  on  the  right  hand  side  can  be  bounded  by 


(A. 19) 


2?i 


|€n(z)I 


|dz|  < 


max  |Cj^(z)|  <  J  max  )Re(  e^(  z )  )  I  ^ 
|z-z|=«  |z-z|-« 


< 


-  Va  , 


max  I Im(e  (z)  )  < 

|z-i|-6 


where  the  last  inequality  follows  from  the  maximum  principle  for  harmonic  functions.  By 
(A. 18-19),  we  therefore  have  |s^f"(z)|  <  Va  .  a  contradiction.  • 
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Proof  of  Thaorm  2.5.  Let  f  e  A(fl) 


For  lone  U  >  0,  to  be  epeclfied  below,  require 


(A. 20) 


0 


-  1 
k-N 


a  q  (E(B)))e^^*de 
k  k 


0,  j  -  M,  M  +  1,  N  t  2, 


Let  bj,  Cj),  be  Fourier  coefflcienta  defined  by 

N  ••  1  f(*{e))e^l*de,  J  -  0,1,2,...  , 

3  211  Q 

1 

‘’jk  2w  (qk(*(e))  -  e*^‘'*)e^l*de,  k  -  0,1,2,...,  J  -  0, ±1,12,...  . 


From  q|j(c(a))  "  *  0(k~^),  and  from  tlie  ■noothnaae  of  qj^tzCs)),  it  followe  that 

Cj^  -  0(k"'’r''),  k  ♦  «,  or  j  ♦  ■  . 

The  linear  ayaten  (A. 20)  haa  the  form 


1  ♦  0(N"^)  0(N"’(>I+1)"’)  0(S“’(N+2)"’) 

“n 

N+1 

0((H+1)"^N“'’)  1+0((N+1)"*)  0(  (H+1)"’(Ht2)"’) 

•n+1 

(A. 21) 

- 

0((H+2)"V’) 

•n+2 

e 

. 

. 

e 

• 

.  ■ 

a 

• 

Ha  obtain 


(A. 22)  bj,+j  -  a„+j  ♦  0((U+J)"’)  I  a^o(k'’>  . 


-  '/y 

By  the  Schwarz  Inequality  the  aum  can  be  bounded  by  0(N  ^  ),  N  -»  ■.  Select  an  N 

aufficlently  large  ao  that  the  linear  ayatem  (A. 21)  la  non-singular.  By  (A. 22) 

*  |b„.4l  ♦  0((N+J)“’),  j  ♦  ».  This  aecurea  the  convergence  of 

1  a^qy(z),  z  e  3n.  Hence,  we  only  need  to  approximate  a  function 
k-H  ^  - 

f*<z)  «“  f(z)  -  I  a.q^(z)  which  ia  analytic  in  0,  and  whose  boundary  values  f*(z(8)) 
K  k 

are  orthogonal  to  e^***,o^  ,e^  ^•**^**, . . ,  .  Similar  to  the  construction  of  f*(z), 


we  can  subtract  suitable  multiples  of  ' 


from  qji  j  “  0(1)N-1,  to 

obtain  q^(z),  j  “  0(1)N-1,  which  all  are  analytic  in  8  and  whose  boundary  values  are 

orthogonal  to  e^^®  ,e^ .  .  Since  the  set  ^®  linearly 

independent,  the  set  ^  1®  linearly  independent.  By  Theorem  2.4, 

f*  e  span  {q*,  q*  ,  ....  q*  ).  This  completes  the  proof  of  the  theorem,  ■ 

o  1  w*1 

Proof  of  Theorem  2.7. 

Let  z(t)  e  [0 ,2ii  1  he  a  parametric  representation  of  88  and  introduce  for 

z  ”  z(s)  e  88 


2t 


P,(zl  :=  J  ln|z  -  z(t)|dt  . 


P^{z)  Is  harmonic  in  (t\(8B  U  {»} ) .  Let  PjCz)  denote  a  harmonic  conjugate  of  F^(z). 
Then 


2« 


F(z)  i“  P^(z)  +  iFjCz)  “  ^  i  ln<*  •  t(t))dt  +  ia  , 


where  a  is  a  real  constant.  Subtracting 


1  2ir 

ln(e^«)  -  —  I  ln(el®  -  e^'^jdt 
0 


gives 


2^ 

F(z(a))  -  ln(e^®)  “  J  ln(-^^-|^ - )dt  +  ia 

2”  i  e^®  -  e**" 


We  choose  a  smooth  branch  of  the  logarithm,  and  by  a  reasoning  analogous  to  that  of 
Theorem  2.1  one  obtains,  with  the  tj  defined  by  (R.5), 


2s 


2s  ,  ,  ,  f  z(B)  -  z(t^) 

z(s)  -  ,  t  ^ 


1  J  ln{ 


e^®  - 


j-1 
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,  It . 

Is  j 

e  -  e 


JL)  .  ’*8 


)(N 


)  , 


N  ♦  <»,  uniformly  for  0  <  s  <  2ii.  This  yields 


exp(N(F(z(s) ) 


ia)  )e 


-INs 


N 

n  (z(b) 
1 


N 

n 

1 


(e 


is 


z(t^)  ) 

It; 


exp(0(N  "  ® ) ) , 


N  »  -  , 


and  by  ( A. 14) , 


n  (z(s)  -  z(tj))  -  exp(N<F(z(3))  -  ia))(1  -  e"^**“)(1  +  0(n"”"®))  . 


Uelng  the  nodes  of  (A.5)  in  the  discretization  gives 

N 

n  (z(s)  -  z(t  ))  -  exp(N(F(z(s))  -  io))(1  +  +  OCn”""®))  , 

1  ^ 


and 


qj,(z(8) ) 


2e~^”*(1  t  0(n'"~^)) 
2(1  +  0(k"""®) ) 


g-iIIS(i  +  oCN"""®)),  N  ♦ 


This  shows  (2.19).  If  z(e)  is  analytic,  we  proceed  analogously  but  utilize  the 
geometric  decay  of  the  discretization  error.  ■ 

Proof  of  Theorem  3.1. 

Conformal  mappings  ‘  'l’j(Wj)  “  ♦jlCg)  *  1  have  for 

z  e  the  representation 

(A. 23)  ln(i|i.(z)(z  -  w.))  -  J  ln(z  -  C )  |S  (z.wj  |dc  |  -  ln(cj  -  ia,,  j  -  1(1)1  , 
j  ^3(1  an  j  )  ) 


where  Glz,v^)  denote  the  Green's  functions  In  (3.1).  The  a^  are  real  constants, 

0  <  Oj  <  2s,  depending  on  the  location  of  Cq.  The  Cj  are  real  constants.  The  value 
of  Cj  is  the  capacity  of  3(5  relative  to  a  charge  at  Wj.  In  (A. 23)  and  below,  we 
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select  smooth  branches  of  the  logarithms.  When  Wj  »  *  the  conformal  mapping  has 


the  representation 


rtf' 

ln(i|i.(2))  “  (  ln(z  -  C)  (C.-KdC)  -  ln(cj  -  la.,  j  -  1(1)1 

i  io  ^  3 


Hie  constant  Cj  in  (A. 24)  Is  the  capacity  of  312.  A  special  case  of  (A. 24)  Is 
(A. 25)  ln(w)  -  J  ln(w  -  C)|dc|,  |w|  >  1  . 

lcl-1 

Average  the  t  equation  (A. 23)  and  use  the  parametric  representation  given  through 
the  density  function  (3.1).  With  ln(c)  ^  In(Cj)  and  ®  “  'J'  ^  ®  j '  obtain 

t  ..  I  1/11^" 

(A. 26)  ln(  n  '  +  ln(  H  ln(  z-z(t)  )dt-ln(c)  -ia,  zeOVatJ, 

j=1  ’  j.1  ^  0 

where  smooth  branches  of  the  tth  roots  are  selected.  For  notational  convenience  we 
assume  that  all  |wj|  <  ..  Let  w  in  (A. 25)  be  defined  by  (3.2),  and  subtract  (A. 25) 
from  (A. 26), 


t  2» 

ln(  n  (Z  -  wj’''*)  j  ln(  ^  -  ln(c)  -  ia 

j»1  ^  ^  0  w(z)  -  e*' 


Discretize  the  integral  in  (A. 27)  by  the  trapezoidal  method  with  nodes  tj  of  (A. 5). 
Lemma  2.1  is  valid  also  for  the  present  parametric  representation  z(t),  and  by  an  error 
estimate  analogous  to  that  in  the  proof  of  Theorem  2.1,  we  obtain  for  z  C  0^130, 

I  /,  1  N  z  -  z(t  ) 

(A. 28)  In  n  (z  -  w.)  “  Z  I  - TT”)  *  O***  "  ln(c)  -  ia,  N  ♦  •>  , 

'  5-’  w(z)  -  e  i 

uniformly  for  z  in  closed  subsets  of  (1.  Multiply  (A. 28)  by  N  and  exponentiate. 


n  (z  -  z(t  ))  ^ 

— - -  n  ltf(2)  -  e  ^)e***'*exp(o(N"  ^  ®)),  N  ♦ 

c"  rt  (z  -  w  i-’ 

j-i  i 
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By  (A. 14)  this  Simplifies  to 


Using  the  nodes  in  the  discretization  of  (A. 27)  yields 


n  (z(s)  -  z(tj) 


N  _  ,  ,N/1 

c  It  (z  -  w .  ) 

1 


(elNs  +  +  o(n"”*®)),  N  ♦  »  , 


n“1 


uniformly  for  0  <  s  <  2it.  This  shows  (3.5).  For  z  e  0^  the  following  representation 
holds 


lt> 


(A. 30) 


ln(  z 


2»  z  -  t 

-  Wj)  =  —  J  ln( - ^ - )dt  -  ln(c  )  -  la 

ft  A  J  -» 


This  can  be  shown  for  example  by  the  calculus  of  residues.  We  write  (A. 30)  as 

_ _  3G 

,  ^  /  c  \  ) 

3(1  e 


(A. 31)  ln(z  -  Wj)  =  J  (CrWj)|dcl  -  In(c^)  -  ia^  . 


Averaging  the  equation  (A. 31)  for  j  *  1(1)1  gives 


(A. 32) 


A 

In  It  (z  -  w  )^^*'  ”  /  ln(^-  )dt  -  ln(c)  -  ia  , 

j  =  1  ^  0  e“ 


where  we  choose  a  smooth  branch  of  the  1th  root.  Discretize  the  Integral  in  (A. 32)  by 
the  trapezoidal  rule  with  nodes  tj.  This  yields 

’’  1/£  1  "  *  "  ***^-v’  -n-2-B 

In  It  (z-w)'  =—  I  ln( - tt: — 2—)  +  0(N  ”)  -  ln(c)  -  ia  , 

j=1  ^  J-1 


as  N  This  convergence  is  uniform  for  z  in  closed  subsets  of  II  .  Multiply  by 

”  it. 

N,  exponentiate  and  utilize  It  e  J  “ 

j”1 


-1.  We  obtain 
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